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We consider tunneling processes in QFT induced by collisions of elementary particles. We propose a semi- 
classical method for estimating the probability of these processes in the limit of very high collision energy. 
As an illustration, we evaluate the maximum probability of induced tunneling between different vacua in a 
(1 + l)-dimensional scalar model with boundary interaction. 



PACS: 11.15KC, 03.70+k, 03.65Sq 

In many models of field theory one encounters tun- 
neling transitions between states separated by an en- 
ergy barrier of finite height Eg, the famous examples are 
false vacuum decay in scalar theories [J and topology- 
changing transitions in gauge-Higgs theory I3 . In 
the weak coupling regime, the rate of tunneling at zero 
energy is exponentially small 0121, but one suspects the 
suppression to vanish once energy exceeding the height 
of the barrier is injected into the system. A particular 
way of inducing the tunneling process is a collision of 
two highly energetic particles. It was conjectured some 
time ago @] that the collision-induced tunneling pro- 
cesses may become unsuppressed at high collision en- 
ergies. Semiclassical study of scalar and gauge-Higgs 
theories |21 E| showed, however, that this is not the 
case: the tunneling probability remains exponentially 
small even if the collision energy E exceeds consider- 
ably the barrier height Es- Furthermore, analyses of 
toy models |3 |H| and unitarity arguments [HI suggest 
that the collision-induced transitions should remain ex- 
ponentially suppressed even if E tends to infinity. In 
other words, it was proposed that, contrary to the ini- 
tial conjecture, the probability of the process has the 
form, 

ViE) oc e-^(^)/s' , (1) 

where is a small coupling constant, and the sup- 
pression exponent F{E) remains positive at all ener- 
gies. This has been confirmed recently by direct cal- 
culation |lUj of the suppression exponent in the whole 
range of energies in a toy two-dimensional model. In 
addition, it was found in |10| that the suppression expo- 
nent reaches its minimum F^ at certain optimal energy 
Eo and remains constant above this energy, 



FiE) = F^, E>E, 



(2) 



One may ask whether or not the formulas (0, (0) are 
valid for other models, with model-dependent values of 
F E 

In this letter we give a general semiclassical method 
for evaluating the minimum value F^^ of the suppression 
exponent and the energy Eq at which this minimum is 
achieved. Our procedure is essentially an adaptation of 
the method of Ref. ^I] , but it is more straightforward 
and technically simpler. We consider the case in which, 
after the appropriate rescaling of the fields, the action 
S takes the form S — S/g'^, where S does not explicitly 
depend on the small coupling constant g. One observes 
that g^ plays effectively the role of the Planck constant, 
and the limit g —^ which we consider below, corre- 
sponds to a semiclassical situation. 

Our starting point is the inclusive probability of tun- 
neling from states with a given number of incoming par- 
ticles and any energy. 



7'.n(iV)=^|(/|Z:^(T/,T,)PjvN)|' 



(3) 



Here U{Tf,Ti) is the evolution operator, and Pn de- 
notes the projector onto states with N particles. The 
initial and final states, \i) and |/), respectively, are at 
the different sides of the potential barrier. Here and 
below the limit Tf +oo, Ti — cx) is assumed. As 
we will shortly see, the quantity ijJJ) can be evaluated 
semiclassically provided the initial number of particles 
is parametrically large, N = N / g"^. We will see that the 
result has a typical exponential form, 



(4) 
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We use the result (0J as a source of information on the 
probability of collision-induced tunneling. It is clear 
that the inclusive multiparticle probability Vni{N) sets 
an upper bound on the two-particle probability of inter- 
est, ViE), at arbitrarily high energies E. Indeed, the 
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energy of the initial state in can be arbitrarily high, 
while any initial two-particle state can be promoted to 
the multiparticle one by adding a number of "specta- 
tor" particles which do not interfere with the tunneling 
process. Hence, the exponential suppression of Vm{N) 
entails the exponential suppression of V{E), and the 
inequality 

Fm(^) < (5) 

holds. Following Rcf. ITj, we conjecture that 



lim Fn,{N) 



Frr 



(6) 



The conjecture ^ is based on the observation that in 
the leading semiclassical approximation the probability 
of tunneling does not depend on the details of the initial 
state provided the initial number of particles is much 
smaller than l/g^. The analogs of the formula © have 
been checked in various situations ^1 El • 

Now, we proceed to the semiclassical evaluation of 
the multiparticle exponent F^{N). Our strategy is to 
represent the inclusive probability Vra{N) in the form 
of path integral and evaluate the latter making use of 
the saddle-point technique. Taking into account that 
Pn ■ Pn — Pn, one recasts Eq. (O in the form, 

VUN)= J2 {m'){^'\PN\^''){^''\m 

V[<j>f,a, a\h, b*]e-^ ^ dKa.al+b^K\ (7) 

X {cl,f\U\a) {a\P^\b) 

where are eigenstates of the field operators which 
we denote collectively by 0, while |a), \h) are the coher- 
ent states, 
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a\Aa) = — la) 



Hereafter we use the shorthand notation U = U{Tf, Ti). 
Making use of the standard path integral for the tran- 
sition amplitude in ^-representation, one writes 

{(Pf\U\a)= j I?(/.exp|l(j5[(/.]-|-i?,(0„a))|, 

<P{Tf) = 4>f 

(8) 

where the boundary term Bi comes from the initial ma- 
trix element ((b\a). 



Bi{(t)i, a) 
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-^0i(k)0j(-k) -t- V2c^kak0j(k)| 



In this expression 0i(k) stands for the spatial Fourier 
transform of the fields at t = Ti. For the matrix ele- 
ment of the projector Pn one obtains (see e.g. |15II16| ) 



{a\PN\b) 



dOexp 



NO 



(9) 

Substituting Eqs. |(HJ), lO into the expression ||7J) and 
performing integration over the variables h, b* , one ob- 
tains the desired path integral representation, 

^m(iV)= j V[cb,cl)',a,a*]dee'^/3\ (10) 

4,(Tf)=<P'{Tf) 

where 



- B,i^,,a) - B*{(j)'^,a) + J dk a^Okc' 



(11) 



Note that the integration over in H1U|) comes from the 
path integral representation for the complex conjugate 
amplitude [{<j)f\U\b)]* . 

The functional F defined in Eq. Hll() is independent 
of the coupling constant g. Hence, at weak coupling 
the integral (|10(l is saturated by its saddle point. The 
saddle-point equations are as follows. Extremization 
with respect to and (f>' gives the classical field equa- 
tions 

SS 6S ^ , 

Boundary conditions for these equations are obtained 
by varying the expression Hll|) with respect to the ini- 
tial and final values of the fields. Using the relation 
5S/6(j)(Tf,x) — (j3{Tf,x.) and taking into account the 
constraint (see. Ea. HlOfl 'l 



one obtains, 



</)(T/,x) = 0'(r/,x) , 



0(r/,x) = (^'(T;,x). 



(12b) 



(12c) 



In the initial asymptotic region Ti —oo the evolution 
of the fields 6, 6' is linear, and one writes^-' 



dk 



(27r)i/2^2^ 



(/kC 



-g_^e M e 



"^'For concrctcncss wc assume that the initial state is an exci- 
tation above the vacuum = 0. 
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^'-J (27r)i/2V2^^ ^^^e + )e . 

The variation of the functional F with respect to 
<j)^, a, a* yields the following relations between the fre- 
quency components, 

/k = /ke^ ^7k=5ke-^ (12d) 

The set of saddle-point equations H12(l can be simplified 
if we recall that the configurations (j) and (j)' saturate the 
amplitude and its complex conjugate, respectively. This 
suggests the Ansatz (f>'{t,x) — [(/)(t,x)]*, which is com- 
patible with the boundary value problem (|12ll provided 
the saddle-point value of 6 is real. Then, the bound- 
ary value problem is formulated in terms of a single set 
of fields x). The conditions (|12b|l . (|12c() imply the 
reality of the fields in the asymptotic future, 

Im(/)(t, x) — > , Im(/)(t, x) — > as t ^ +oo , (13a) 

while Eqs. Ijl2d(l read 

/k = e- V . (13b) 

The boundary condition in the asymptotic past, 
Eq. I)13b(l . can be understood as follows. In the limit 
6 — !■ +00 it coincides with the Feynman boundary 
condition and thus corresponds to the initial state with 
semiclassically small number of particles, 0; finite 

9 picks up the most favourable state with non-zero N. 

The number of equations in the boundary value 
problem (|12a|l . p3(l is equal to the number of unknowns. 
Generically, for a given value of 6 this problem has a 
unique solution (prtit,^)- We call this solution "real- 
time instanton", as it lives on the real time axis, in 
contrast to the ordinary instanton which is defined in 
Euclidean time. Note that the boundary condition 
(|13b(l implies that the real-time instanton is complex- 
valued. On the other hand, its imaginary part should 
vanish in the asymptotic future due to the condition 
()13a|l . Let us discuss the consequences of this prop- 
erty. Assume that at large (but finite) times the solu- 
tion gets linearized about some real static configuration, 
(f>rt{t, x) — (f>s{x) + S(j){t, x). Then Eq.(12a) implies that 
Im50 — > as t — !■ -|-oo. This amounts to requiring that 
the configuration ips is unstable, so that ImSfj) evolves 
along its negative mode, ImScj) oc e~'^°"''' *. The natural 
candidate for (ps (x) is the static solution "sitting" on top 
of the potential barrier separating the sectors of initial 
and final states. Accepting the terminology of gauge- 
Higgs theories 0| we call this solution "sphaleron"'^^ 

■^'in case of scalar theories such solution is known as "critical 
bubble" Q], 



We arrive at the conclusion that the real-time instan- 
tons describe formation of the sphaleron as < — > +oo. 
This is in fact a common property of solutions relevant 
for collision-induced tunneling at energies higher than 
Es ^1 El El ■ The transition is completed by the de- 
cay of the sphaleron into the states of interest, which 
proceeds with probability of order one. 

The last saddle-point equation obtained by varying 
the functional pi|l with respect to the parameter 9, re- 
lates the value of 9 to the initial number of particles, 

I dkf^g^ . (14) 

Substituting the solution into Eq. Hll|l . one obtains the 
formula with the suppression exponent 

F^{N) = ~N9 + 2 Im I^S[(brt] + \j dx^rfi^n | ■ 

'""^'(15) 

Note that the term in braces is the action of the real- 
time instanton integrated by parts with respect to time. 

At finite value of the parameter 9 (non-zero value 
of N) the real-time instanton is, by construction, a 
smooth solution to the equations of motion and has a 
well-defined classical energy Eo{N). At t ^ +00 the 
real-time instanton describes semiclassical evolution in 
the final state. Thus, Eo{N) coincides with the energy 
of the final state saturating the probability (PJ. This 
means that if we restrict the sum in (jJJ to the states 
with the fixed energy E = Eo{N) (cf. Ref. illj) we ob- 
tain the same result (0}. One concludes that Eo{N) is 
the optimal energy for tunneling from the states with 
given number of particles N — N / . The limit 

E, = hm E,{N) (16) 

determines the optimal energy for collision-induced tun- 
neling^^. It is easy to understand the most favourable 
transition at higher energies. The system releases the 
energy excess {E — Eq) by a perturbative emission of a 
few particles (which costs only a power suppression in 
g^) and tunnels at the optimal energy Eq [3E|. Due 
to this process F{E) stays constant at _E > E^, i.e. the 
formula Q holds. 

Let us illustrate our method by considering a simple 
example. We consider free massless scalar field (j){t, x) 
living in (1+1) dimensions on a half-line a; > 0, with 



^'in general, we cannot exclude the situation when Ea(N) — > 00 
as A' — > 0. In the example below, Eo is finite. 
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self-interaction localized at the boundary point a; = 0. 
The action of the model is^^ 



S 




dx{df,4>y - - cos((/)(t, 0))] 



(17) 

where the parameter /z sets the characteristic energy 
scale of the boundary interaction. The model p7|l is 
used in solid state physics to describe the transport 
in quantum wires ^^'^ Josephson chains with de- 
fects ^]. A detailed semiclassical treatment of the 
model (|17|l is given in 10 . 

The model (|17|) has a number of vacua 0„ — 27m, 
n = 0,±1,..., which are separated by the potential 
barriers of the height Es = determined by the 

maximum of the boundary potential. One also finds 
an unstable static solution (jj^ = tt, sphaleron, which 
"sits" on the top of the "first" potential barrier. The 
process we are interested in is tunneling between the 
vacua = and (j) = 2tt induced by a highly energetic 
particle scattering off the boundary. To calculate the 
minimum suppression exponent and the optimal energy 
of this collision-induced tunneling process, one finds the 
family of real-time instantons by solving the boundary 
value problem (|12a|l . (|13|l . 

Since the bulk evolution of (j){t, x) is that of free 
massless scalar field, we represent the general solution 
in the form 

(l){t,x) = ipin{t + x) + ipont{t - X) , (18) 

where ipin and i^Sout are the incoming and outgoing 
waves. The boundary interaction leads to non-linear 
equation at a; = 0, 



dx(l> = /isin0, X — 



(19) 



Due to the condition H13a|l . the outgoing wave ipout is 
real. Introducing the real and imaginary parts of the in- 
coming wave, finiO = o(?)+*^(C)j one rewrites Eq. lfTn|) 
as a set of two real equations, 



b' = /i sh 5 cos u , 

u' = 2a' — /i ch 6 sin u 



(20) 
(21) 



where u(^) = a(f) -|- 'Pout(0- The remaining conditions 
in the asymptotic past, Eq. (|13b|l . should be imposed 



on the frequency components of the incoming wave ipi, 
To this end, one performs Fourier expansion of ipi^, 



ipUt + x)= / dfc ^i„(fc)e^'=(*+-) , 



(22) 



and finds that the positive and negative frequency com- 
ponents of the solution, and g'Lf., are proportional 
to (pia{—k) and ipin{k), fc > 0, respectively. Thus, 
Eq. Ijl3b|l takes the form fin{—k) = e"^ [ipin{k)]* , where 
fc > 0. It is straigforward to check using the Cauchy for- 
mula that the latter condition is equivalent in its turn 
to the following relation between the real and imaginary 
parts of the initial wave (see Ref. jlUj). 



1 



1 



1 

— V.P. 



da^^ , (23) 



where the integral is understood in the sense of princi- 
pal value. We use Eq. H23() as an alternative formulation 
of the condition (|13b|l in our model. 

One also requires that the functions ipin, (font have 
appropriate asymptotics. To ensure the finite number 
of particles in the initial state we require the incom- 
ing wave to be well localized, ifmiO — > as ^ — > ±00. 
Besides, as the initial state is an excitation above the 
vacuum (j) — 0, we write (/?out(C) — >■ at f — > —00. On 
the other hand, as we have already discussed, the rele- 
vant solution contains the sphaleron at t ^ +00. Thus, 

</'out(C) -> TT as ^ -> +CO. 

The problem (|20|) . H21|) . H23|l can be solved numeri- 
cally by the following iterative method. At each cycle 
of iterations one starts from the function^) u — u(°)(^) 
and solves Eq. (^0)) explicitly. 



biO = Inth |^-|^%osu(")(ei)d6 



(24) 



where >r is an integration constant. Then, integrating 
numerically Eqs. (|^ . if^ . one finds a(^) and m(^) for 
given value of x. Finally, one picks up the value of x 
such that the function u{^) has correct asymptotics as 
^ — > ±00. In this way the improved approximation for 
u{^) is obtained, and a new cycle of iterations begins. 
After 30 cycles one obtains the solution with precision 
of order 10"^. 

Given the family of real-time instantons, one calcu- 
lates numerically the suppression exponent F^(0) via 
Eq. I|15|) . and also the energy Eo{9). To calculate the 



®' Normally, one should use some infrared regularization for the 

massless scalar model in (1-1-1) dimensions. However, the specifics ^'At the very first cycle one takes, e.g., u^^^(S^) = n/2 + 

of regularization turn out to be irrellevant for our purposes. arctan(/i5). 
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FIG. 1. Suppression exponent (left) and energy of real- 
time instantons (right) versus the number of incoming 
particles in the model (1171 . 

number of particles it is convenient to recast Eq. H14|) 
in the form 

N^g^N^^^ J d^a'iOm ■ (25) 

Functions Fn^{9), Eo{9) and N{9) determine the depen- 
dence of the suppression exponent and the energy on 
N. These are plotted in Fig.^ Taking the Hmit iV ^ 
one obtains Fm = 10.27, Eq — 2.&hEs- Thus, the semi- 
classical suppression factor is exp{—lQ.21/g^) at all en- 
ergies exceeding Eo = 2.65Es, and the suppression is 
even stronger at lower energies. 

In conclusion we summarize our method. To cal- 
culate the minimum suppression exponent F^ and the 
optimal energy Eq of a collision-induced tunneling pro- 
cess, one finds the family of complex classical solu- 
tions, real-time instantons, satisfying the boundary con- 
ditions H13() . Given the real-time instantons, one calcu- 
lates the number of particles N, the suppression expo- 
nent Fm(iV), Eqs. (Hlj), pHll, and the energy EoiN). By 
itself, the quantity Fi^(N) provides the lower bound on 
Era, Eq. lO; the Hmit N ^ yields Fm and Eo accord- 
ing to Eqs. ®, lfTH|l . 
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